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ABSTRACT 

We analyze the high- resolution X-ray spectrum of Hercules X-l, an intermediate-mass X-ray binary, 
which was observed with the XMM- Newton Reflection Grating Spectrometer. We measure the elemental 
abundance ratios by use of spectral models, and we detect material processed through the CNO-cycle. 
The CNO abundances, and in particular the ratio N/O > 4.0 times solar, provide stringent constraints 
on the evolution of the binary system. The low and short-on flux states of Her X-l exhibit narrow line 
emission from C VI, N VI, N VII, O VII, O VIII, Ne IX, and Ne X ions. The spectra show signatures of 
photoionization. We measure the electron temperature, quantify photoexcitation in the Hea lines, and 
set limits on the location and density of the gas. The recombination lines may originate in the accretion 
disk atmosphere and corona, or on the X-ray illuminated face of the mass donor (HZ Her) . The spectral 
variation over the course of the 35 d period provides additional evidence for the precession of the disk. 
During the main-on state, the narrow line emission is absent, but we detect excesses of emission at 
~10-15 A, and also near the O VII intercombination line wavelength. 

Subject headings: X-rays: binaries — line: formation — line: identification — pulsars: individual (Her 
X-l) — accretion, accretion disks — binaries: eclipsing 



1. INTRODUCTION 

Hercules X-l is a bright intermediate-mass X-ray binary 
which has been observed extensively after its discovery by 
Schreier et al. (1972). The system contains an X-ray pul- 
sar with Ppuisc = 1-24 s period (Tananbaum et al. 1972; 
Wilson, Scott, & Finger 1997). A synchrotron resonance 
feature yields a magnetic field of B = 3.5 x 10 12 G (dal 
Fiume et al. 1998; Truemper et al. 1978). Some mod- 
els yield a lower magnetic field (Baushev & Bisnovatyi- 
Kogan 1999). Optical light-curves (Gerend & Boynton 
1976) and X-ray eclipses (Tananbaum et al. 1972) yield a 
Porb = 1-7 day orbital period. The 1.5 ± 0.3 M Q neutron 
star has a 2.3±0.3 M & companion, HZ Her, which changes 
from A to B spectral-type over the orbital period, due to 
the strong X-ray illumination on its surface (Reynolds et 
al. 1997). Pre-eclipse and anomalous dips in the X-ray flux 
are observed (Giacconi et al. 1973; Scott & Leahy 1999) 
due to the interaction of an accretion stream with the ac- 
cretion disk. The unabsorbed luminosity of Her X-l is 
L = 3.8 x 10 37 ergs -1 , using a distance of D = 6.6±0.4 kpc 
(Reynolds et al. 1997). 

Hercules X-l exhibits an unusual long-term X-ray flux 
modulation with P* = 35 d period (Giacconi et al. 1973). 
Changes associated with this period have also been ob- 
served in the optical light curves (Gerend & Boynton 
1976), X-ray pulse shapes (Deeter et al. 1991), X-ray dips 
(Giacconi et al. 1973; Scott & Leahy 1999), and X-ray 
spectra (Ramsay et al. 2002). The 35 d X-ray light-curve 
is asymmetric and contains two maxima: a state of ~ 8 d 
duration reaching the peak flux F max named the main-on, 
and a secondary high state of ~ 4 d duration reaching 



~ 1/3 -Fmax named the short-on. A low-flux state with 
~ 1/20 F m ax ensues at other epochs. The period P* varies 
from cycle to cycle, and it has been observed to be 19.5, 
20, 20.5, and 21 times P or b in the course of > 5 yr of 
continuous monitoring with the All-Sky Monitor onboard 
the Rossi X-ray Timing Explorer (Scott & Leahy 1999). 
The "J = phase is defined as the time the main-on state 
begins. Observations indicate that "J" = coincides only 
with orbital phases = 0.23 or 0.68 (Scott & Leahy 1999). 

The 35 d cycle has been associated with a tilted ac- 
cretion disk that precesses by some unknown mechanism 
(Petterson 1975; Petterson, Rothschild, & Gruber 1991; 
Scott, Leahy, & Wilson 2000). The X-ray light-curve, the 
variations in the pulse profiles, and the variability of the 
dips with ^P-phase are fit by a geometric model of a pre- 
cessing, warped accretion disk, together with the beam of 
a pulsar (Scott, Leahy, & Wilson 2000). Among the pa- 
rameters obtained in this fit are an 85° inclination with 
respect to the line of sight, a 20° precession opening an- 
gle for the outermost disk, and an 11° precession angle 
for the innermost disk. The latter precession angles are 
dependent on assumed input parameters such as the disk 
thickness. 

The ultraviolet (UV) spectrum exhibits line emission 
from species such as C IV, N V, and O V, which orig- 
inate in two separate components that produce superim- 
posed broad and narrow lines. The intensity of the narrow 
lines has the same orbital variation as the UV continuum, 
which is thought to originate on the illuminated face of 
HZ Her. The broad line component follows roughly the 
velocities expected for the accretion disk (Boroson et al. 
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1997) . Observations during eclipse ingress and egress sug- 
gest that the broad line region originates in a prograde 
disk of - 10 11 cm radius (Boroson et al. 2000). The UV 
lines (i.e. N V) have a weak P Cygni profile component 
which indicates the presence of a wind outflow (Chiang 
2001). 

The Her X-l broadband X-ray emission can be described 
as a blackbody component with temperature kT ~ 90 eV, 
plus a power-law component with a 24 keV exponential 
cutoff, and a 42 keV cyclotron feature (dal Fiume et al. 

1998) . An Fe K fluorescence line evolves with \I>-phase. 
The XMM-Newton EPIC data, presented by Ramsay et 
al. (2002), show a 6.4 keV Fe K line which is practically 
unresolved during the low and short-on states, and a broad 
line at 6.5 keV with 330 ±20 eV FWHM during the main- 
on. 

In this article, we present the high-resolution spectrum 
of Her X-l in the 5 to 38 A band, which we obtained 
with the XMM-Newton Reflection Grating Spectrometer 
(RGS). The high-resolution X-ray spectrum evolves dra- 
matically with 35 d phase. We analyze the spectral data 
from three distinct epochs (section 2). The spectrum dur- 
ing the low and short-on states is a faint power-law con- 
tinuum plus many narrow emission line features from a 
photoionized plasma (section 3). In contrast, the main-on 
state has a bright continuum, neutral absorption features, 
excess of emission near O VII Flea, and a continuum which 
cannot be fit well with blackbody and/or power-law mod- 
els (section 4). The observed emission lines arise from the 
recombination of electrons with ions and from the subse- 
quent radiative cascades. The gas is predominantly heated 
by photons due to the large luminosity and small volume 
of the system. 

We analyze and model the high- resolution spectra of Her 
X-l to obtain the elemental abundance ratios and measure 
the state variables of the gas, and we make preliminary 
identifications of the newly detected emission region(s). 
We apply spectral diagnostics on the plasma to constrain 
its location (section 3.1), its density (section 3.2), and to 
measure the electron temperature (section 3.3). We intro- 
duce a method to extract the elemental abundance ratios 
from the low and short-on spectra. The method involves 
extracting the emission measure distribution from model 
fits of the recombination line fluxes (section 5) . We discuss 
whether the narrow line emission originates in a disk at- 
mosphere and corona, on the illuminated companion, or in 
both (section 6.1). We describe the implications of the de- 
tection of CNO-processed material for the evolution of the 
X-ray binary (section 6.2). We outline a physical scenario 
which produces broad emission lines, which is associated 
with the magnetopausc and the inner disk in X-ray pulsars 
(section 6.3). In that context, we speculate on the origin 
of the complex continuum observed during the main-on 
(section 6.4). Our conclusions are in section 7. 

2. OBSERVATIONS 

Three observations with XMM-Newton RGS during the 
main-on, short-on and low states, reveal X-ray line emis- 
sion that varies with <F-phase. The RXTE All-Sky Moni- 
tor (ASM) public data are used to find the epoch of main- 
on start, or ^ = (see the light curve in Fig. 1). The 
values of the cjy and ^-phases during the XMM-Newton 



observations are shown in Table 1. 

We describe the processing of each RGS spectrum. All 
the data are processed with the most recent version of the 
science analysis software (SAS v. 5.3). Response matri- 
ces are extracted using 4000 energy bins in the full RGS 
energy range. Hot pixels which are flagged by the pipeline 
processing are excluded from the analysis. Continuum 
and line fits are performed with the XSPEC v. 11.1 anal- 
ysis software (Arnaud 1996). The total RGS background 
count-rate is 0.09 c s _1 in the low state, 0.08 c s _1 in 
the short-on state, and < 0.28 c s _1 in the main-on state. 
The background is roughly independent of wavelength, at 
a level of <~ 2.7 x 10~ 3 c s _1 A -1 in the low and short- 
on states. The background level rises by a factor of 2 to 
4 for wavelengths A < 7 A. The background level in the 
main-on is very small (< 0.015 c s _1 A -1 ) compared to 
scattering in the cross-dispersion direction from first or- 
der photons. The background-subtracted event lists from 
the SAS pipeline are used throughout. The extraction re- 
gions from the RGS data are inspected and left unchanged 
from the SAS default. During main-on, RGS 2 is set to 
high timing resolution mode, as part of the performance 
verification of the instrument, while RGS 1 is set to spec- 
troscopy mode. The RGS guaranteed time observations 
were performed during the low and short-on states, with 
both RGS in spectroscopy mode. We present the spec- 
troscopy mode data. Since CCD chip 7 in RGS 1 and CCD 
chip 4 in RGS 2 are not operational, the wavelength ranges 
of 10.3 to 13.6 A and 19.9 to 24.9 A are excluded from each 
RGS dataset, respectively. For the main-on data, we have 
enough counts to use the RGS 1 second order, which fills 
the wavelength gap in the first order. 

3. SPECTROSCOPY OF THE LOW AND SHORT-ON STATES 

The low-state RGS spectrum (shown in Fig. 2) con- 
tains line emission from Nc IX, O VIII, O VII, N VII, 
N VI and C VI ions. Within the resolution, no veloc- 
ity broadening has been observed. The continuum shape 
is described by a power-law, although the goodness-of-fit 
is likely degraded by a collection of weak and unresolved 
spectral lines. There is a neutral O edge, which is mostly 
instrumental. The line and continuum measurements are 
shown in Tables 2 and 3, respectively. No evidence for 
absorption local to the source is observed on these states, 
implying that the line emission region is optically thin to 
both lines and continuum. 

The total RGS count-rate in the short-on is two times 
larger than during the low-state (also in Fig. 2). There 
is variability within each observation as well, but we leave 
that out of our analysis. The line luminosities are usually 
larger in the short-on than in the low state (Table 2). The 
better statistics during the short-on allows us to detect the 
O VII and N VII radiative recombination continua (RRC). 
On Figure 3, we show the added spectrum of the low and 
short-on states, which have similar features, to clarify the 
statistical significance of the detection of Ne IX Hea, C VI 
Lya, the O VII RRC, and the N VII RRC (the line pro- 
files in this figure are not adequate for analysis because 
we have added spectral data with different spectral reso- 
lution). The EPIC 0.3-10 keV count rates, which are 398, 
54 and 21 c s _1 in the main-on, short-on, and low states, 
respectively (Ramsay et al. 2002), indicate fluxes which 
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are consistent with the long term XTE ASM light curve 
(Scott & Leahy 1999), with our main-on observation near 
peak, and our short-on observation at half-peak. 

Our phenomenological model fit to the spectrum con- 
sists of a power-law continuum, emission lines, absorption 
edges, and RRC. To perform the continuum fits, and to 
use the x 2 statistic, we bin the spectrum with a minimum 
of 20 counts per bin. The spectra shown in Figures 2, 
3, and 4 are rcbinncd as indicated on each caption. The 
spectral lines are fit at full resolution with Gaussian func- 
tions by independently varying their normalization, width, 
and centroid. The line identifications are shown in Table 
2. The continuum is re-fit for accuracy after the line fits 
are added. The ISIS analysis software (Houck & Denicola 
2000) is used to check the line fits from XSPEC. With 
ISIS, the continuum level is measured in two regions ad- 
jacent to the lines. The line equivalent widths (EW) cal- 
culated from that continuum level are shown in the last 
two columns of Table 2. The line fluxes from XSPEC and 
the equivalent widths measured with ISIS agreed to < 10 
%, although the difference is - 25 % for N VII Hea, N VI 
Hea, and O VIII Lya. The difference is due to some hot 
pixels which are removed only with XSPEC, and to the 
uncertainty of the continuum flux level. 



3.1. Helium-like ion emission line diagnostics 

Electron density and temperature measurements in the 
10 8 < n e < 10 18 cnr 3 and 10 6 < T < 10 7 K ranges 
can be obtained with helium-like ion emission line ratios 
(Gabriel & Jordan 1969; Porquet & Dubau 2000). The 
dominant ionization mechanism can be determined from 
the same lines (Liedahl 1999). Low-Mass X-ray binaries 
(LMXB), and other X-ray binaries with accretion disks, 
should exhibit densities high enough to render these diag- 
nostics useful. 

In Her X-l, the helium- like ion line complexes show a 
prominent intercombination line blend (i denotes the two 
transitions 2 3 Pi, 2 — > l 1 ^) and a weaker resonance line (r 
denotes 2 1 P\ — > l 1 ^), as shown in Figures 5 and 6. The 
forbidden line (/ denotes 2 3 Si — > l 1 5o) is not detected 
for any of the helium-like ions Ne IX, O VII, and N VI. 
The notation n 2S+1 Lj identifies each two-electron atomic 
state, with quantum numbers n, L, S, and J. 

The line ratio G = (f+i)/r is consistent with ~ 4 for all 
cases shown in Table 4, which indicates that photoioniza- 
tion is the dominant heating mechanism. For collisionally 
ionized gases, G < 1, decreasing for larger temperatures 
(Liedahl 1999; Porquet & Dubau 2000). The G ratio in- 
dicates the electron temperature is T < 2 x 10 6 K, or 
kT < 170 eV (Porquet & Dubau 2000). 

We set upper limits for the R = f/i line ratio (see Tables 
2 and 4). The R ~ ratio indicates that 1) the density is 
above the threshold calculated by Porquet & Dubau (2000) 
and Bautista & Kallman (2000), thereby exciting the 2 3 Si 
level to the 2 3 Pi j2 levels by electron impact, and/or 2) the 
2 3 Si level is photo-excited to the same 2 3 Pi ;2 levels, as de- 
scribed by Mewe & Schrijver (1978), and later identified 
in an O-star wind by Kahn et al. (2001). The density lim- 
its implied by the measured R, ignoring photoexcitation, 
would be n e > 4 x 10 11 cm~ 3 for N VI, n e > 2 x 10 12 cm~ 3 
for O VII, and n e > 5 x 10 12 cm~ 3 for Ne IX (Porquet & 



Dubau 2000). However, the photoexcitation rate is 

2 

ire 

Wf->i = F v f osc , (1) 

m e c 

where e, m e are the electron charge and mass, c is the 
speed of light, / osc is the oscillator strength, and F v is 
the flux (in photons s _1 cm~ 2 Hz -1 ) at frequencies res- 
onant with the 2 3 Si — ► 2 3 Pi i 2 transitions, in this case in 
the UV. We use the / osc calculated by Cann & Thakkar 
(1992). The calculated rates (Table 4) imply that UV pho- 
toexcitation can compete with electron impact excitation, 
such that the R ratio density diagnostics are not directly 
applicable in Her X-l. 

Instead, we use the R ratio to set limits on the dis- 
tance between the X-ray line region and the UV source, d. 
Because the radiation field is diluted with increasing dis- 
tance, Wf_»i oc F v . i oc 1/g? 2 . Beyond a threshold distance 
d cr it, which is set by u>f~>; = uif, the photoexcitation rate 
Wf^i is smaller than the radiative decay rate uif, so the 
/ line is not affected by photoexcitation, and neither is 
R. For d < d CT it, photoexcitation can suppress the / line 
and enhance the i line, until the R ratio approaches zero 
when the UV radiation field is sufficiently strong. To cal- 
culate rf cr it, we use the distance to Her X-l of D = 6.6 kpc 
(Reynolds et al. 1997), the UV flux F v , t measured by 
Anderson et al. (1994) and Boroson et al. (1996), and we 
assume a UV point source. Our results are shown in Ta- 
ble 4, which yield d clit ~ 10 12 cm. Kahn et al. (2001) 
used this sort of procedure for a collisionally ionized gas. 
The value of G? cr ;t becomes ~ 1/2 of that in Table 4 if we 
use the F{\ = 1637 A) - 0.4 x 10~ 13 erg s" 1 cmT 2 A" 1 
measured by Boroson et al. (1997) at (j> = 0.66. This 
flux is in rough agreement with observations by Boroson 
et al. (2000) at <p = 0.685-0.764. We attribute this to 
time- variability associated with both <fi- and 'J'-phase, and 
possibly luminosity fluctuations. 

Boroson et al. (2000) found that the UV flux peaked 
at orbital phase <p — 0.5, and interpreted this as the il- 
luminated face of HZ Her. We can use this information 
to set a limit to the distance between the pulsar and the 
narrow line region (rii nc ). Since rii ne < \J <i 2 rit + rf A , and 
cfcrit > fhi ~3x 10 11 cm, we estimate rij no < d cr ;t, where 
is the distance between Her X-l and the LI Lagrange 
point. This limit on rii no , obtained from photoexcitation, 
is shown by the vertical line in Figure 7. 

3.2. Limits on the density from thermal and ionization 

balance 

Because the R ratio is insufficient to quantify the den- 
sity, we use the ionization parameter £ instead. The ion- 
ization parameter is defined as £ = L/n e r 2 erg cm s _1 
(Krolik, McKee, & Tarter 1981), where L is the luminosity, 
n e the electron density, and r is the distance from the neu- 
tron star to the X-ray line emitting region. The ionization 
state of an optically thin, photoionization-dominated gas 
in equilibrium is fully determined by £, and by the spec- 
tral shape of the ionizing continuum. Assuming thermal 
and ionization balance, and with the main-on continuum 
as input, we calculate the charge state distribution for a 
grid of £ by use of the XSTAR plasma code (Kallman & 
McCray 1982). Our XSTAR models and recombination 
emission calculations (section 5) determine the £ where 
the recombination flux from a given ion is significant. For 



4 



Jimencz-Garate et al. 



a fixed emission measure (section 6.1), the £- values which 
correspond to the FWHM of the recombination line power 
of Nc IX, O VII, and N VI Hea lines are log 10 £ =1.3-2.0, 
1.0-1.9, and 0.9-1.7, respectively. 

To partially account for the existence of optically thick 
gas in Her X-l, we do not make density estimations with 
£ and instead set limits which are more robust. We ac- 
count for the decrease in number of photons available for 
photoionization, for the case where 1) the line emission 
region is optically thin, and 2) the line-of-sight from the 
line emission region to the neutron star may be optically 
thick. We do not account for line emission in an opti- 
cally thick region, hence we do not need a true transfer 
computation. One large source of uncertainty in the den- 
sity estimates obtained from £ is our knowledge of the flux 
impinging on the line emission region. In Her X-l, the 
X-ray flux variation with ^-phase is associated with large 
changes of the optical depth in our line-of-sight to the pul- 
sar. Since the optical depth between the pulsar and the 
line emission region can be r > 0, we set upper limits to 
the density from n e — Le~ r /£r 2 < L m ax/£minf 2 , where 
L max is the main-on state, or maximum, X-ray luminos- 
ity. We take £ m i n = 1.0 for O VII from the T lower limit 
set in section 3.3, and from the £(T) solution from XS- 
TAR. Because £ m i n coincides with the lower end of the 
FWHM of the O VII Hea recombination line power, our 
calculations are validated by the data. We use £ m ; n and 
L max = 3.8 x 10 37 (L>/6.6 kpc) 2 erg s" 1 to obtain the den- 
sity upper limits for O VII shown in Figure 7. The density 
limits for Ne IX and N VI are ~ 0.50 and ^1.3 times 
those for O VII, respectively. As the density decreases be- 
low the limit, r increases and the £ of the emission region 
increases. 

3.3. Radiative Recombination Continuum 

Free electrons recombining with ions produce an RRC 
width proportional to T. The observed RRC photons cor- 
respond to recombination to the K shell. For a photoion- 
ized gas, the RRC will be narrow, since T -C /, where I is 
the ionization energy (Liedahl & Paerels 1996). 

The N VII RRC at A = 18.59 A is detected with a > 10 a 
confidence, yielding (6.5 ± 2.8) x 10 -5 photons cm~ 2 s _1 , 
and kT = 3.7 ± 1.3 eV (90% confidence). We also detected 
the O VII RRC at A = 16.78 A with a - 4 a confidence 
level in the combined short-on and low state data. The 
O VII RRC is visible in the added spectrum of the low 
and short-on states (Fig. 3). The O VII RRC normal- 
ization is (2.8l2'J) x 10~ 5 photons cm~ 2 s" 1 , and kT is 
2-6 eV (this 90% confidence interval is used to set a density 
limit in section 3.2). The RRC temperatures confirm pho- 
toionization is the dominant heating mechanism (Liedahl 
& Paerels 1996). Other RRC are not detected. 

The observed N VII RRC is not far from the expected 
flux, if compared to N VII Lya. The N VII ion has the best 
measured RRC-to-Lya ratio, which is 0.28 ± 0.14. The 
hydrogen-like RRC should have ~ 0.75 the flux (Liedahl 
& Paerels 1996) of the corresponding Lya flux (see Table 
2). Thus, the N VII RRC is close to, but weaker, than the 
expected strength. For both O VII and N VII, the mea- 
sured T agree with our XSTAR calculations, within the 
statistical errors. 

A high-tcmperature component, or a finite optical 



depth, can explain the weak RRC observed in the narrow 
line region. An RRC with T > 50 eV, which may originate 
in an extended hot corona (section 6.1), would be difficult 
to distinguish from the continuum. Alternatively, if the 
optical depth is significant, RRC photons produce ioniza- 
tions and "fill in" the absorption edges, thereby pumping 
the recombination lines. However, optical depth effects 
also imply changes in the Lya/Ly/3 ratios, and perhaps on 
the Hea r and Lya lines, which are not observed. 

4. SPECTROSCOPY OF THE MAIN-ON STATE 

The main-on RGS spectrum is characterized by a bright 
continuum which has a broad "bump" at ~10-15 A (Fig. 
4). Excess emission near the O VII Hea wavelength is 
present. We performed a single power-law fit on the spec- 
trum, using the first order data, and we filled the wave- 
length gap with the second order data. 

The RGS spectrum is dominated by a soft continuum 
component which is distinct from the hard (photon index 
= 0.87) power-law component observed with EPIC (Ram- 
say et al. 2002). This hard power-law component does not 
significantly contribute to the RGS spectrum for A > 8 A. 
For simplicity, we fit the RGS continuum level indepen- 
dently (Table 3). The EPIC PN is consistent with the 
RGS spectrum at 5 < A < 14 A, within the known uncer- 
tainties in the calibration. 

Centered at the O VII i line wavelength, we observe an 
excess of emission (see Fig. 5). Due to the noise in this 
feature, we are not certain if it is composed of a single line 
or multiple lines, or if it is an emission line at all. How- 
ever, rcbinning reveals a statistically significant emission 
feature. For reference, we perform a two-Gaussian line 
fit, assuming the feature is O VII Hea, with R — and 
G = 4 (the line ratios observed in the other states). This 
yields a a = 3200 ± 800 km s _1 velocity broadening, but 
we caution that a may not have a physical meaning. We 
search for a feature at the O VIII Lya wavelength, but the 
statistics are not sufficient for a detection (see the observed 
profile in Fig. 8). We set upper limits for the O VIII Lya 
and N VII Lya line fluxes, assuming they have the same 
velocity broadening as the putative O VII Hea. 

We detect a neutral O edge and set a limit on the neutral 
N edge (see edge data in Table 5). The RGS calibration 
data in the SAS version 5.3 analysis software does include 
an instrumental O (oxide) edge. Cold absorption models 
with solar abundances do not provide acceptable fits. We 
fit the O edge to obtain an upper limit on the interstellar 
O I column density of No i < 4.8 x 10 17 cm~ 2 , which im- 
plies an equivalent H column of (Na) cq < 5.7 x 10 20 cm~ 2 , 
if we use the Wilms, Allen, & McCray (2000) solar abun- 
dances. A resonant ls-2p absorption feature from atomic 
O is detected. This ls-2p interstellar absorption feature 
has been observed before in other X-ray binary spectra 
(Paerels et al. 2001). We set a 90% upper limit for the 
column A N i < 1.2 x 10 17 cm~ 2 . The Ne and Mg edges 
are not detected, but limits cannot be set reliably due to 
the uncertainty in the continuum in the ~ 10-15 A range. 

Analysis of the first and second order RGS 1 spectrum 
does not yield any robust line identifications of individ- 
ual Fe L transitions. A power-law fit results in negative 
residuals at 9-10 A, and positive residuals at 11-15 A and 
~ 7-8 A(Fig. 4). These residuals do not appear thermal, 
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and they can be described as a broad "bump" at ~10- 
15 A. The "bump" at 10-15 A was previously identified 
with Fe L emission (Oosterbroek, Parmar, Martin, & Lam- 
mers 1997), but these observations were performed with 
low-resolution instruments. The high- resolution spectrum 
obtained by the Einstein spectrometer also showed a broad 
emission feature at 1 keV with low statistical significance 
(McCray et al. 1982; Vrtilek et al. 1991). The possible 
nature of this "bump" is discussed in section 6.4. 

5. A METHOD TO MEASURE ELEMENTAL ABUNDANCE 
RATIOS FROM RECOMBINATION EMISSION LINES 

The Hea and Lya emission lines observed in the low and 
short-on states are modeled to extract the elemental abun- 
dance ratios between C, N, O, and Ne. We calculate the 
C VI Lya, N VI Hea, N VII Lya, O VII Hea, O VIII Lya, 
and Nc IX Hea line fluxes from a grid of regions which have 
£ separated by A(log 10 £) = 0.1. Good ^-sampling is nec- 
essary for an accurate abundance determination because 
the recombination rates and charge state distribution are 
sensitive to £. To calculate the state of the gas in each 
region, we run the XSTAR code, using as input the main- 
on broadband continuum shape. The main-on continuum 
consists of a 92 eV blackbody, a broken power law which 
has a break at 18 keV and photon indeces 0.88 and 1.8, 
and a 24 keV exponential cutoff. This spectrum was mea- 
sured by dal Fiume et al. (1998), and it is consistent with 
the EPIC observations. We assume thermal and ionization 
balance, which is justified due to the short time required 
for equilibration at the densities of interest (see the limits 
on the density in section 3.2). Knowledge of the emission 
measure EM — J n\dV as a function of £ is necessary to 
measure the elemental abundances, which implies a model 
dependence. To constrain EM(£) from the spectrum, one 
may choose a single £ for each ion, such that the radia- 
tive recombination flux is maximized, and then calculate 
the EM for all the ions detected (Sako, Liedahl, Kahn, & 
Paerels 1999). However, each ion species emits copiously 
at a broad range of £, which must be taken into account 
to obtain an accurate measure of EM{^) and the abun- 
dances. 

We introduce a method to find the emission measure 
and the elemental abundances simultaneously. We assume 
a functional form for EM(£) and then test the likelihood of 
the hypothesis with \ 2 - The differential emission measure 
distribution is well fit with the power-law form 

The zero-slope d(EM) /d(log 10 £) = 0, and the commonly 
assumed d(EM)/d£ = K <5(£ — £„) (this is equivalent to a 
single ionization parameter), are inconsistent with the ob- 
served spectrum, where K, 7, and £ are fit parameters and 
6() is the delta function. We use the recombination rate 
coefficients calculated by Liedahl (private communication) 
with the HULL AC code (Klapisch et al. 1977). The ob- 
served fluxes are fitted using a simplex \ 2 -minimization al- 
gorithm, with 5 free parameters: the relative abundances, 
K, and 7. We fit each observation independently, to test 
the robustness of the method. The results in Table 6 show 
that 1) consistent fits can be found by the algorithm, 2) 
the phenomenological equation (2) provides an excellent fit 
to the RGS spectra in the low state and a marginal fit in 



the short-on state, 3) the results are sensitive to the shape 
of the ionizing spectrum, which can be seen by comparing 
the results obtained with the main-on spectrum and with 
a 20 keV bremsstrahlung spectrum, and 4) the N/O abun- 
dance ratio is at least four times the solar value, while the 
C/O and Ne/O ratios indicate that both carbon and oxy- 
gen are depleted with respect to neon. The marginal fit of 
the short-on spectrum may be due to a physical emission 
measure which differs from a power-law distribution. The 
large (N VII + N VI) /(O VIII + O VII) ratio is evident in 
the spectrum on Figure 2. 

We discuss the validity and limitations of the method. 
The method is robust insofar the DEM(^) is well-covered 
in £ by the observed ions, and because the emission mea- 
sure is not extrapolated to other £. Also, it relies on the 
observed absence of optical depth or resonance scatter- 
ing effects in the helium- like ion line ratios. The abun- 
dance ratios depend on the XSTAR plasma code assump- 
tions, which may introduce systematic errors that are 
not addressed in this paper. We have used a 20 keV 
bremsstrahlung as a benchmark to test the dependence 
of our results to the ionizing spectrum, and we do not 
observe an improvement in the x 2 of the fit; albeit the 
C/O and N/O abundance ratios are lower than with the 
main-on spectrum. XSTAR is run at n e — 10 10 cm -3 , 
a regime for which three-body recombination of the ob- 
served ions is unimportant, and in the optically thin case. 
For n e < n cr it, the recombination line fluxes depend only 
weakly on n e . The value of n cr ;t depends on the ion, and is 
roughly n cr i t ~ 10 17 cm~ 3 for O VIII (Kallman & Bautista 
2001). The low-density assumption holds for the narrow 
line region (Fig. 7). We assume Gaussian statistics since 
there are > 44 photons in each of the lines used for the 
abundance ratio measurement. We excluded Ne X Lya 
from the fit since it is detected with low statistical signifi- 
cance, but our fits are consistent with the observed limits. 

Models of the Her X-l UV spectrum roughly agree with 
the abundances we measure. CNO-processing abundances 
were also evidenced by the N V/O VI line ratio observed 
in the UV (Boroson et al. 1997), which was compared to 
spectral models of a photoionized accretion disk (Raymond 
1993). We crudely estimate an abundance ratio N/O ~ 6.0 
times solar from the observed discrepancy of the N V/O VI 
line flux ratio with the CNO-disk model. However, the 
N V/O V ratio yields N/O ~ 2.6 times solar. A C IV/O V 
ratio depleted below the solar abundance value was also 
observed by Boroson et al. (1997). Our X-ray spectral 
measurements are more reliable than the UV ones because 
they are less sensitive to radiation transfer effects, and be- 
cause the X-ray line fluxes have been fitted consistently 
with the DEM. 

The functional form of DEM(t), through the K and 
7 parameters, contains information on the gas density 
and geometry. A power-law dependence of the DEM(£) 
is expected in some simple geometries such as: a uni- 
formly filled sphere, an optically-thin disk atmosphere 
with a scale height H oc r n , and an illuminated slab with 
p(z) oc p e~ z / z °. A full interpretation of the K and 7 
parameters is outside the scope of this paper. 
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6. DISCUSSION 

A new, X-ray narrow line component in the Her X-l sys- 
tem has been discovered. We discuss the possibility that 
this narrow line region is the illuminated face of HZ Her, or 
an accretion disk atmosphere and corona. Our elemental 
abundance measurements serve to quantify the CNO pro- 
cessing products from HZ Her. We discuss the connection 
between the CNO products and the history of the binary 
system. 

The main-on spectrum shows a peculiar continuum fea- 
ture and possibly a broad emission line component. We 
discuss the physical plausibility of such a broad line com- 
ponent. 

6.1. Recombination emission from the disk atmosphere 
or the illuminated companion 

We investigate the nature of the photoionized gas emit- 
ting narrow X-ray lines during the low and short-on states. 
The UV line emission consists on a narrow (cr ~ 80 km s _1 ) 
line component, which is likely to originate on the illumi- 
nated face of HZ Her, plus a broader (cr ~ 300 km s _1 ) 
component, attributable to the accretion disk (Boroson et 
al. 2000). Similar dynamical components in the X-ray nar- 
row line region would not be distinguished with the RGS 
energy resolution alone. 

The limits found on the density and location of the nar- 
row line region are shown in Figure 7. We set the density 
lower limit of the O VII region with the emission measure 
calculated from the line flux. The lower limit to the root- 
mean-square of the density is (n rms ) m ; n = y 'EM /F max , 
where the maximum volume U max — ^7rr 3 . The i?-ratio 
density thresholds listed in section 3.1 indicate that, within 
much of the allowed region in Figure 7, the 2 3 S — > 2 3 P 
electron-impact excitation rate is larger than the photoex- 
citation rate. Radius limits for the putative Kepler orbits 
are obtained from the circularization radius and from up- 
per limits to the velocity broadening (Fig. 7). However, a 
disk larger than the circularization radius (~ 1.7x 10 11 cm) 
was deduced from UV and optical light curves (Howarth 
& Wilson 1983). 

An extended, photoionized accretion disk atmosphere 
and corona may be responsible for the observed recom- 
bination emission and the underlying continuum. Model 
calculations (Jimenez-Garate, Raymond, Liedahl, & Hai- 
ley 2001) show that a centrally illuminated, photoion- 
ized accretion disk would develop an extended, Compton- 
temperature coronal structure, and a more compact, X- 
ray recombination-emitting atmospheric layer (see Fig. 9). 
We first interpret the behavior of the Her X-l spectrum in 
terms of this picture. The line fluxes observed are within a 
factor of two (scaled by the Her X-l luminosity and disk ra- 
dius) of the photoionized accretion disk atmosphere model 
results, and the line ratios are all within a factor of unity 
of model values. The model density of the helium-like ion 
region is 10 13 -10 14 cm -3 if the disk radius is r ~ 5 x 10 10 
cm, which are inside the limits on Figure 7. If the X-ray 
line velocity broadening can be resolved, this would sup- 
port the accretion disk emission hypothesis (the evidence 
is marginal from the RGS, see Table 2). The line emit- 
ting disk atmosphere is optically thin towards most lines 
of sight, although its flattened geometry allows it to par- 
tially shield itself from the pulsar (e~ T < 0.1). The soft 



X-ray continuum is dominated by Compton scattering in 
the disk corona during the low and short-on states. The 
optical depth of both atmosphere and corona decreases as 
the 1) radius increases, and as 2) the height above the 
disk mid-plane increases. At high inclination, the inte- 
rior of the disk is obscured by the outer disk edge. As 
the disk precesses out of its edge-on configuration, parts 
of the disk which are close to the neutron star are less 
obstructed, and deeper layers of the disk come into view, 
increasing both the atmospheric line flux and the Comp- 
ton scattered light from the corona. This may explain 
the spectral evolution observed from the low to the short- 
on state. As the disk inclination decreases further in the 
main-on state, the pulsar comes into view and overwhelms 
the narrow line flux from the disk. The duration of X-ray 
eclipse ingress or egress at the low and short-on states is 
measured in the 1500 s to 3 hr range (Leahy 1995; Choi, 
Scon, Dotani, & Nagase 1997), corresponding to source 
radii of 4 x 10 10 cm to 3 x 10 11 cm. Eclipses of the X-ray 
lines arc very likely to occur, since similar radii estimates 
have been obtained from the eclipse of the broad UV line 
emission region. About ~ 0.1 x 10 -13 erg s _1 cm~ 2 A -1 or 
~ 6 % of the UV continuum, together with the broad C IV 
and N V emission lines, are gradually eclipsed, indicating 
a ~ 10 11 cm UV emission region (Boroson et al. 2000). 
Eclipse measurements, together with the X-ray spectrum 
and its variability, and the photoionized disk models, lend 
support to the disk atmosphere and corona origin of the 
X-ray emission during the low and short-on states. 

The spectrum of the accretion disk corona (ADC) source 
4U 1822-371 (Cottam et al. 2001) is very similar to the Her 
X-l low and short-on state spectra, confirming that ADC 
sources provide an edge-on view of the disk and supporting 
the idea that the accretion disk in Her X-l is precessing. 

The variability of the X-ray line emission with ^-phasc 
and 0-phase, combined with the predicted disk inclina- 
tion from the Scott, Leahy, & Wilson (2000) model, sug- 
gests that at least some of the X-ray line flux originates 
on HZ Her. The total line emission flux in the RGS 
band at phase ^ = 0.60 is two times larger than dur- 
ing phase $ = 0.26, at approximately the same orbital 
phase (4> <~ 0.5). Taking the predicted outer disk inclina- 
tion from the Scott, Leahy, & Wilson (2000) model, we 
have drawn the schematic in Figure 10 (Model B). If the 
line emission originates in an accretion disk atmosphere 
and corona, the line flux should increase as more of the 
disk's area comes into view. If the line emission originates 
in HZ Her, the line flux should increase as the visible area 
of HZ Her increases. If the outer disk inclination predicted 
by the model of Scott, Leahy, & Wilson (2000) is correct, 
then some of the line emission in the short-on must be 
coming from the photoionized face of HZ Her. If the X- 
ray line emission originates on HZ Her, it should vary with 
orbital phase in a way similar to the UV continuum flux 
(Boroson ct al. 2000). Some dependence of the narrow 
line flux on (/>-phase is hinted by the main-on observation. 
The upper limits to the narrow lines which we derive for 
the main-on observation at <j> <~ 0.2, are below the line 
fluxes in the low and short-on states, which were observed 
at <f> ~ 0.5. However, this may be due to a ^-phase de- 
pendence only. Observations at various orbital phases and 
similar ^P-phase, during the low and short-on states, can 
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disentangle the companion and disk components. 

However, model values of the outer disk precession phase 
differ (see Models A and B in Fig. 10). Model B (Scott, 
Leahy & Wilson 2000) yields a precession phase which lags 
by AW ~ 0.16 the phase from model A (Gerend & Boynton 
1976; Howarth & Wilson 1983), which fits the optical and 
UV lightcurves. In model A, the maximum disk opening 
angle to our line-of-sight occurs during the main-on and 
short-on peaks, at <!/ = 0.13 and 0.63. Model A favors the 
interpretation that the line emission originates in a disk 
atmosphere and corona, since the projected disk opening 
angle and the line fluxes during the short-on are larger 
than during the low state. 

6.2. New constraints on the X-ray binary evolution from 
the derived CNO abundances 

The unusually high N/O abundance ratio and the low 
C/O and O/Ne ratios we have measured are the result of 
H-burning by the CNO-cycle in the core of a massive star. 
The CNO-cycle produces N and depletes O and C, until H 
is exhausted (Clayton 1983, and references therein). Evi- 
dence for CNO processing in the Her X-l system was previ- 
ously obtained from UV spectra, which show N-enrichmcnt 
and C-depletion (Boroson et al. 1997, see section 5). Evo- 
lution models of an isolated star with a mass comparable 
to HZ Her cannot explain the observed abundances, even 
after the end of the H-burning phase. A mechanism to 
enrich the HZ Her envelope with CNO-processed material 
is required, and such a mechanism is likely linked to the 
presence of a companion. The CNO abundances thus pro- 
vide a stringent constraint on the evolution of the binary 
system. 

We compare the measured CNO yields with those in an 
isolated intermediate-mass star. Stellar evolution models 
(Schaller, Schaerer, Meynet, & Maeder 1992) of an isolated 

2.5 Mq star show a N/O abundance ratio enhancement of 

2.6 to 2.8 times the zero-age N/O ratio, a C/O ratio of 
0.60 to 0.70 times its zero-age value, and an unchanged 
Ne/O ratio. The quoted range corresponds to zero-age so- 
lar (Z = 0.02) and sub-solar (Z — 0.001) metallicities, at 
the end of the H-burning phase, after ~ 5.8 x 10 8 yr. The 
N-yield and the C-depletion are maximized near the end 
of the H-burning phase. The N-yield increases with stellar 
mass. For example, a 4.0 M Q star produces a N/O ratio 
which is 4.0 times its zero-age value (Schaller, Schaerer, 
Meynet, & Maeder 1992). Models of the evolution of an 
isolated ~ 2.5 M star under-predict the observed N-yield 
in HZ Her. 

The observed N-enrichment with RGS cannot be at- 
tributed to systematic uncertainties in the data, nor could 
we attribute it to a dependence on the photoionization 
equilibrium calculations. The systematic errors due to fit- 
ting of the emission lines and the continuum are < 25% 
(section 3). Any RGS calibration errors are ~ 10 % where 
most of the emission lines have been measured, and even 
the worst-case O I edge calibration error scenario implies 
a < 35% error in the O VII Hea line flux, since to i ~ 0.3. 
If we propagate a ±35% uncertainty in the O VII Kca flux 
through the abundance extraction calculation, the N/O 
ratio changes by < 35%. The measured abundance ratios 
depend on the shape of the ionizing spectrum, which af- 
fects the equilibrium state of the plasma (section 5). The 



changes in the 0.5-10 keV continuum from the main-on 
state to the low and short-on states are well-fit with a 
varying Nh on a constant-slope power-law (Ramsay et 
al. 2002). Most of the line emission should be generated 
by reprocessing of the direct pulsar (main-on) continuum, 
which consists of a power-law with a cutoff. Taking a 
20 keV bremsstrahlung ionizing spectrum still yields a 
N/O > 4.0 times the solar value, and underabundant C 
and O with respect to Ne (Table 6), so a systematic dif- 
ference in the plasma equilibrium state cannot, by this 
measure, explain the CNO abundances. 

The detection of CNO-processed material is validated 
by measurements of the (C + N + 0)/Ne abundance ra- 
tio. The (C + N + O) abundance sum is conserved by the 
CNO-cycle, so the sum should be near the solar value. We 
measure the (C + N + 0)/Ne ratio for the two best fits in 
Table 6, obtaining 0.74±0.28 times solar with the main-on 
continuum model, and 0.45 ± 0.19 times solar with the 20 
keV bremsstrahlung model. 

To explain the observed abundances in HZ Her, a mecha- 
nism is required to transfer CNO-processed material to the 
envelope. One possibility is that enhanced mixing and sig- 
nificant mass loss from the HZ Her envelope enrich it with 
CNO products from the core. The stellar core is richer in 
CNO products: for the 2.5 M Q example above, the abun- 
dance ratios are C/O - 0.04, N/O - 10, and Nc/O - 10 
times the zero-age values (Schaller, Schaerer, Meynet, & 
Maeder 1992). These abundances have not reached the 
equilibrium CNO values, which have higher N-yields. The 
CNO products must have been efficiently transported from 
the core to the surface. However, the stellar core contains 
only ~ 0.1 times the mass of the star. The shedding of a 
N-poor stellar envelope though Roche-lobe overflow may 
explain the observed N-enrichmcnt. Intermediate-mass 
companions may survive super-Eddington accretion, as ev- 
idenced by observations of the Cygnus X-2 mass donor 
(King & Ritter 1999), and by models of the evolutionary 
track of that binary system (Podsiadlowski & Rappaport 
2000). The mixing of the core's CNO products with the en- 
velope gas is needed to explain the observations, but that 
may complicate the calculation of the CNO burning rates, 
which depend on the concentration of nuclear species in 
the core. Another possibility is that the CNO-processed 
material was transferred from the neutron star progenitor 
onto HZ Her. A massive neutron star progenitor can read- 
ily produce the CNO elements in a ~ 10 7 yr timescale, 
and would not require HZ Her to be highly evolved. The 
presence of a companion has likely affected the elemental 
composition of HZ Her, which presents a challenge to our 
understanding of the evolution of the binary. 

It is highly likely that the evolutionary state of the 
companion stars in low- and intermediate-mass X-ray bi- 
naries is more evolved than had previously been antici- 
pated. This has been shown by modeling the evolution 
paths of 100 sample systems, from a grid of initial stellar 
masses and periods (Podsiadlowski, Rappaport, & Pfahl 
2002). These models show that many LMXB evolve from 
intermediate-mass binaries, so that many LMXB should 
be H-deficient and He-enriched, and the surface compo- 
sition of the evolved secondaries should show evidence of 
CNO processing. The CNO abundances we find in Her 
X-l are consistent with this basic picture. The determi- 
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nation of the age of HZ Her, its zero-age mass, and the 
initial period of the binary can be better addressed by bi- 
nary evolution models now that the CNO abundances have 
been measured (Podsiadlowski, Rappaport, & Pfahl 2002, 
and private communication). 

6.3. The physical plausibility of a broad line component 

We make the case that a broad line region is physically 
plausible if it is neighboring the Her X-l magnetopause, 
and we point out that broad recombination lines have been 
observed previously in another accreting X-ray pulsar. 

A low-mass, compact X-ray binary system, 4U 1626-673, 
contains a pulsar which is accreting from a white dwarf, 
and it exhibits prominent recombination emission lines 
(Schulz ct al. 2001). The Chandra High-Energy Transmis- 
sion Grating detected the Ne X Lya, Ne IX Hea, O VIII 
Lya, and O VII Hea lines with FWHM ~ 2500 km s _1 , 
and it resolved them into two distinct peaks, suggesting a 
disk origin. If the broadening is due to the orbital velocity, 
the material would be located on the disk, exterior to the 
pulsar magnetopause (Schulz et al. 2001). 4U 1626-673 
has a cyclotron line with 0.85 times the energy than the 
corresponding line in Her X-l, a 3.4 times larger corotation 
radius, and broad lines with ~ l/20th the flux of Her X-l 
(the unabsorbed X-ray flux from 4U 1626-673 is <~ l/36th 
the Her X-l main-on flux). It is possible for broad recom- 
bination lines to be associated with the accretion flows 
around X-ray pulsars, which have B ~ 10 n -10 12 G. 

In Her X-l, from the putative O VII Hea broad lines, 
we estimate that the emission region would be inside the 
4.5 x 10 8 cm < r < 3 x 10 9 cm boundary We assume the 
velocity broadening is due to circular Kepler orbits, and 
that ~90% of the emission in the O VII Hea Gaussian line 
fit is enclosed within v = ±2cr of the centroid. The esti- 
mated inner disk radius is just larger than the corotation 
radius r co = (GMP p 2 ulsc /47r 2 ) 1/3 = 1-9 x 10 s cm, and is 
near the Alfven radius (Ghosh & Lamb 1979). 

The putative O VII Hea broad line luminosity would 
imply that a fraction of order ~ 10 _1 of the luminous en- 
ergy of the pulsar is directed towards the line emission 
region. Emission measure calculations allow us to set a 
n e > 10 14 cm~ 3 density lower limit, from the radius upper 
limit above. This implies that there would be significant 
line opacity effects, which would have to be included in the 
photoionization equilibrium calculations. Hypothctically, 
a nearly Compton-thick region enshrouding the pulsar, of 
radius of a few times 10 8 cm, may be illuminating an inner 
accretion disk, which then emits the broad recombination 
lines (Fig. 9). To illustrate, a shell of Compton-thick 
gas at r ~ 10 s cm, and with Ar ~ 10 7 cm thickness re- 
quires n e ~ 10 17 cm~ 3 and log 10 £ ~ 4, which indicates 
the plasma would be fully ionized, except perhaps for iron. 
An illuminated disk region of radius r~3x 10 9 cm would 
have n e ~4x 10 17 cm~ 3 , which corresponds to the £ for 
maximum line power, and taking Ar ~ 0.5r and volume 
V = 2nrhAr, its height would be h ~ 10 3 cm. 

The duration of eclipse ingress and egress varies dur- 
ing the Her X-l main-on state, and it yields source size 
estimates in the 5 x 10 s cm to 3 x 10 9 cm range, albeit 
the smaller estimate may be affected by the HZ Her atmo- 
sphere (Leahy 1995; Leahy & Yoshida 1995; Day, Tennant, 
& Fabian 1988). 



6.4. The 10-15 A "bump" 

The origin of the 10-15 A "bump" observed in the main- 
on state is unclear. We describe three possible explana- 
tions for this "bump". The first possibility is that the 
"bump" is a superposition of Fe L emission lines with a 
velocity broadening as in O VII Hea. The blended tran- 
sitions may range from Fe XVII through Fe XXIV, at 10- 
17 A. However, photoionized gases have weak Fe L emis- 
sion relative to hydrogen-like and helium-like ions of low- 
and mid-Z elements (Liedahl et al. 2001). If the "bump" 
is Fe L emission, it could be due to shock-heated gas in 
the magnetopause, although the G line ratio from O VII 
Hea is not consistent with collisional heating. 

The EPIC data suggest that the "bump" results from 
the overlap of the pulsar power-law continuum with a re- 
processed, soft component. Shortward of A ~ 6 A, for 
the hard power-law, the pulse profile consists of sharply 
peaked pulses, indicating a < 10 s cm emission region. 
Longward of A ~ 14 A, the pulse profile is nearly sinu- 
soidal, indicating a ~ 10 9 cm emission region (Ramsay et 
al. 2002). The soft component is not a Planck spectrum. 

A third possibility is that we are observing a blended 
and broad Fe L resonance absorption complex. Fe L ab- 
sorption is not suppressed in photoionized gases. In this 
case, absorption lines from hydrogenic and helium- like ions 
should also be observed. 

7. CONCLUSIONS 

We analyzed the high-resolution X-ray spectra of Her 
X-l obtained with the XMM-Newton RGS from three ob- 
servations performed at distinct states of its 35 d cycle. 

We detect narrow recombination emission lines during 
the low and short-on states of Her X-l. Emission lines are 
detected from C VI, N VI, N VII, O VII, O VIII, Ne IX, and 
perhaps Ne X, plus weak RRC of O VII and N VII, which 
indicate the gas is photoionized. The velocity broaden- 
ing of the lines is at the resolution limit of the RGS (e.g. 
a < 260 km s" 1 for N VII Hea). In the 5-38 A band, the 
low and short-on states have power-law continua, while the 
main-on continuum exhibits excess emission in the ~ 10- 
15 A band. The continuum flux during the short-on is 
twice as during the low state, a tendency that is followed 
by most line fluxes as well. 

We measure the abundance ratios among C, N, O, and 
Ne with the following method: we use a power-law emis- 
sion measure distribution, in conjunction with a grid of 
XSTAR plasma models, and the HULLAC recombina- 
tion rates, to fit the narrow line fluxes, and obtain the 
abundance ratios with the emission measure parameters 
simultaneously. We perform the measurements for the low 
and short-on state data separately, and we obtain consis- 
tent results with acceptable \ 2 (Table 2). The enrichment 
of nitrogen relative to oxygen, of more than four times 
the solar values, plus the depletion of C and O with re- 
spect to Ne, indicate extensive H-burning in a massive 
star. The measured abundances require a mechanism for 
transferring the CNO-processed material onto the HZ Her 
envelope. This mechanism is likely linked to the presence 
of a companion, providing an additional constraint on the 
evolution of the X-ray binary. 

We use spectroscopic analysis and models to set lim- 
its on the density and location of the narrow line region. 
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These limits (Fig. 7) provide important clues on the na- 
ture of the line emission region. We assume thermal and 
ionization balance to set upper limits on the density, and 
we use the emission measure derived from the line fluxes 
to set lower limits on the density. If the line velocity 
broadening is due to Kepler motion, we may set bounds 
to the orbital radii. We use the Hea line ratios and UV 
photoexcitation calculations to set upper limits to the ra- 
dius enclosing the region. The low RRC temperatures 
(30, 000 < T < 60, 000 K) allow us to validate our density 
upper limits and the photoionization equilibrium models. 

The narrow line region may be identified with an ac- 
cretion disk atmosphere and corona, or with the illumi- 
nated face of HZ Her. The evidence for the disk identi- 
fication relies on the modeled structure and spectra from 
a photoionized disk (Jimenez-Garate, Raymond, Liedahl, 
& Hailey 2001), which agree with the limits set on the 
density (10 13 -10 14 cm -3 ). The unresolved velocity broad- 
ening indicates the outermost (r ~ 5x 10 10 cm) radii of the 
disk dominate the emission, in agreement with the mod- 
els, while fluxes of the observed hydrogen and helium-like 
lines arc within a factor of two of the calculations. Obser- 
vations from previous missions of the eclipse ingress and 
egress reveal an emission region which can match the size 
of the disk and its corona. The flux variation with orbital 
phase and with 35 d phase, on the other hand, favors a 
contribution from the illuminated face of HZ Her to the 



narrow line emission. 

The variability of the Her X-l spectrum lends support to 
the precession of the accretion disk. Regions with contrast- 
ing dynamical properties are coming into view at different 
35 d phases. Notably, Her X-l exhibits an ADC-like spec- 
trum during low and short-on states, indicating an edge-on 
disk, while during main-on states, the spectrum is domi- 
nated by the continuum, due to a disk inclination which 
exposes the pulsar to our line of sight. 

We detect an excess of emission centered on the O VII 
Hea i line wavelength during the main-on state. If this fea- 
ture is due to the O VII Hea i and r lines, their velocity 
broadening would be in the order of 10 3 km s _1 . Similarly 
broad, yet double-peaked lines, have to date only been ob- 
served from the accreting pulsar in 4U 1626-67 (Schulz et 
al. 2001). The observed emission emission feature is noisy, 
but if it is real, it may be due to emission in the inner 
accretion disk, near the pulsar magnetopause. 
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Fig. 1. — Light curve obtained with the RXTE all-sky monitor. 
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Fig. 2. — High-resolution spectra obtained with RGS 1 (red) and RGS 2 (black), and phenomenological fits (thick-lines). The low state 
(13 ks) and short-on state (12 ks) observations are shown on adjacent panels. The spectra have been rebinned adaptively. 
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14 15 16 17 18 19 20 

Wavelength (A) 




Wavelength (A) 

Fig. 3. — Summed RGS spectra of both low-state and short-on data, which highlight the observed RRC of O VII and N VII, and C VI 
Lya. The bin size is 21 and 27 mA for the top and bottom frame, respectively. The line spread function is different for each RGS and is not 
corrected on these spectra. 
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Fig. 4. — High-resolution spectrum with RGS 1 for the main-on state observation (11 ks), plus a phenomenological model fit (thick-line). 
Both first and second order spectra are shown in the upper and lower part of the figure, respectively. 
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Wavelength (A) Wavelength (A) Wavelength (A) 

Fig. 5. — The evolution of the O VII emission line complex with 35 d phase. The intcrcombination (i) and resonance (r) lines are present 
in the low and short-on states. A complex, broad structure during the main-on is plotted at full resolution. RGS 1 data (symbols) and fits 
(solid lines) are shown. 
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Fig. 6. — The evolution of the N VI emission line complex with 35 d phase. The intercombination (i) and resonance (r) lines are present 
in both low and short-on states. RGS 1 data (black symbols), RGS 2 data (red symbols), and the respective fits (solid lines), are shown. 
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Radius enclosing VII narrow line region (cm) 

Fig. 7. — Confidence limits (90%) on the density and position of the narrow line emission region, obtained from spectroscopy and modeling. 
The limits shown correspond to O VII, and similar limits with N VI and Ne IX can be set. 
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wavelength (A) 

Fig. 8. — Main-on observation near the O VIII Lyo wavelength of 18.97 A, at full resolution. Further observations are needed to determine 
whether there is a P Cygni feature. 
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mil 

10 cm 



Fig. 9. — Hypothetical sketch of an accreting, magnetized neutron star, loosely based on the Ghosh & Lamb (1979) model. The accretion 
disk structure is also sketched, based on models by Jimenez-Garate, Raymond, Liedahl, & Hailey (2001). We define the atmosphere as a thin 
skin on the disk which is emitting X-ray recombination emission. The disk atmosphere and corona, or the illuminated face of the companion 
(not shown), may be responsible for the narrow line emission. Plasma from the inner disk can be channeled by the magnetic field, forming a 
highly ionized and possibly Compton-thick structure in the magnetosphere. This structure and the inner disk are potential sites for Compton 
scattering, X-ray absorption, and recombination emission. The tilt of the pulsar magnetic moment with respect to its spin axis is not shown. 



20 



Jimcncz-Garate et al. 




Model A Model B 



Fig. 10. — Schematic of the Her X-l geometry during our three observations, showing the outer disk precession angle according to Model 
A (Gerend & Boynton 1976; Howarth & Wilson 1983), and to Model B (Scott, Leahy, & Wilson 2000). The diagram docs not include disk 
warp, is roughly to scale, and is meant to show the projection of the system to our line of sight. 
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Table 1 

Observations with the Reflection Grating Spectrometer 



MJD J 


XMM 


O bscrvation 


Exposure 


Orbital 


35 d 


RGB 1 


RGS 2 


Interval 


Orbit 


Start 


(ks) 


Phasc (1) 


Phase b 


(c/«) c 


(c/sr 


51935.036-51935.169 


207 


2001 Jan 26. 00:51:13 UT 


11.31 


0.18-0.26 


0.17 (main-on) 


11.29 




51972.927-51973.083 


226 


2001 Mar 4, 22:14:18 UT 


13.41 


0.47-0.56 


0.26 (low state) 


0.37 


0.42 


51984.906-51985.053 


232 


2001 Mar 16, 21:45:01 UT 


12.63 


0.52-0.60 


0.60 (short-on) 


0.83 


0.93 



"Mean Julian Day. 

^Thc zero of the 35 d phase is obtained at MJD 51929.08 (cf> olb = 0.68) and MJD 51964.00 (cf> olh = 0.22) from RXTE ASM dwell-by-dwcll 
light curve data. The ASM light curve shows some evidence (a 4<r detection in one dwell) that the second turn-on occurred at MJD 51963.08 
('/'orb = 0.68) instead, but is affected by pre-eclipse dips. 

c First order count rates. 

References. — (1) Deeter et al. 1998. 



Table 2 

Observed X-ray emission lines 



Linc(s) 


A (A) 


A (A) 


<t (km/s) 


Flux a 


State 


EW (eV) b 


Counts c 


C VI Lya 


33.74 


33.72 ± 0.01 


< 380 


6 4+ 2 ' 1 


Low 


2.7 ± 0.4 


47 






33.75 ± 0.02 


< 320 


10 + 3 


Short-on 


1.8 ± 0.3 


66 


N VI Hcq 


29.53 (/) 






< 0.94 


Low 


< 0.6 










< 1.6 


Short-on 


< 0.4 






29.08 (i) 


29.07 ± 0.01 


270 ± 100 


29 ± 4 


Low 


18.6 ± 1.9 d 


173 d 






29.10 ± 0.015 


< 260 


40 ± 6.5° 
< 30 


Short-on 
Main-on 


10.9 ± 0.9 d 


217 d 




28.78 (r) 


28.77 ± 0.01 


< 260 


6.6 ±2.0 


Low 








28.75 ± 0.02 


< 260 


9 ±5° 


Short-on 






N VII Lya 


24.78 


24 77+0-02 

1 1 -0.01 


< 590 


15 ± 3 


Low 


5.2 ± 0.7 


58 






24.80 ± 0.02 


520+ 230 


23 ± 5.5° 


Short-on 


4.8 ± 0.5 


111 






24.85 ± 0.17 


3200 (fixed) 


84 ±30 


Main-on 






O VII Hcq 


22.10 (/) 




< 3.6 


Low 


< 2.8 










< 3.6 


Short-on 








21.80 (i) 


21.79 ± 0.01 


< 400 


31 ± 5 


Low 


23.4 ± 2.4 d 


154 d 






21.82 ± 0.02 


3200 ± 800 


70 ± 10 
155 ± 30 


Short-on 
Main-on 


22.4 ± 1.6 d 


287 d 




21.60 (r) 


21.61 ± 0.03 


< 1600 


10 ± 4 


Low 








21.62 ± 0.02 


3200 ± 800 


8 ± 4 
38 ± 8 


Short-on 
Main-on 






O VIII Lya 


18.97 


18.96 ± 0.01 


390 ± 200 


12.0 ± 2.5 


Low 


7.5 ± 0.8 


98 






18.96t°;°? 


< 450 


14 ±3 


Short-on 


3.0 ± 0.3 


85 






19.05 ± 0.05 


P Cygni? 


< 90 


Main-on 






Ne IX Hcq 


13.70 (/) 






< 1.6 


Low 












< 6.2 


Short-on 








13.55 (i) 






5-6i 2 2 ;5 


Low 
Short-on 


10.1 ± 1.7 d 
7.2 ± 1.0 d 


44 d 
61 d 




13.45 (r) 






< 4.8 


Low 














Short-on 






No X Lya 


12.13 


12.20 ± 0.04 




7 ± 3 
3.6 ± 3.0 


Low 
Short-on 







Note. — Errors are 90 % confidence limits unless otherwise noted. Symbols: A — theoretical wavelength; 
A — measured wavelength; a — standard deviation; State — Flux state within the Her X-l 35 d cycle; EW = 
equivalent width. Note A, a and the flux arc obtained from Gaussian fits. 

a In units of 10 -5 phot cm _2 s -1 . The errors do not include the continuum level uncertainty, which can be 
significant for the weakest lines. 

b Quoted errors in this column arc 1-er (68 % confidence limits). 

c After pipeline processing and application of response matrices, these are the corrected X-ray counts per line 
(or line complex) with the local continuum subtracted. The EW is then calculated. 

d These measurements include both intercombination (i) and resonance (r) lines. 

c Wc use only RGS 2 data in these line fits, to avoid a few flagged pixels in RGS 1. The RGS 1 still detects 
these lines at levels comparable to RGS 2. 
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Table 3 

Continuum power-law fit parameters 



State 


Photon Index 


Normalization at 1 keV 


N H 


X7DOF 






(phot keV" 1 cm" 2 s" 1 ) 


(10 18 cm 2 ) 


Low 


1.65 ±0.06 


(4.1 ±0.1) x 10" a 


< 86 


562/467 a = 1.21 


Short-on 


1.75 ± 0.03 


(9.8 ± 0.2) x 10~ 3 


< 75 


1179/945 a = 1.25 


Main-on 


1.88 ±0.10 


0.130 ±0.007 


< 57 


4312/3122 = 1.39 



Note. — Errors are 90% confidence limits. DOF = degrees of freedom. 
"The DOF reflect rebinning performed for plotting and fitting the continuum. 



Table 4 

Helium-like line diagnostics and photoexcitation of the 2 3 5' level 





Her X-l 


R = f/i 


G=(j + i)/r 




Flux 


Wf 


iL> f ^ i (r=10 11 cm) 


dcrit 


Ion 


State 


Line 


Line 


(A)' 1 ' 


F X{ ^. (10- 13 erg 


2 3 S -> 1 X S 


2 3 S -> 2 3 P 


Radius 






Ratio 


Ratio 




cm" 2 s- 1 A" 1 ) 


Rate (s- 1 )^ 


Rate (s" 1 ) 


(cm) 


Nc IX 


low 


< 0.3 


> 1.3 


1270 


1.9 ± 0.1 (3) 


1.09 X 10 4 


5 x 10 5 


7 X 10 11 




short-on 


< 0.5 


o o + 16 
d.8_ 2 ! 












O VII 


low 


< 0.14 


3.1 ± 1.4 


1637 


1.5 ± 0.2 (3) 


1.04 x 10 3 


1 x 10 6 


3 x 10 12 




short-on 


< 0.06 


8.7 ± 4.5 












N VI 


low 


< 0.04 


4.4 ± 1.5 


1906 


1.5±0.1 (4) 


2.53 x 10 2 


2 x 10 6 


9 x 10 12 




short-on 


< 0.05 


4.4 ± 2.6 













Note. — Errors are 90% confidence limits. Symbols: Af^; = wavelengths of 2 3 S —* 2 3 P transitions; «Jf_»i = radiative decay rates; u)f_>i = 
calculated photoexcitation rate. 

References. — (1) Porquet et al. 2001; (2) Drake 1971; (3) from the Hubble GHRS at <f> = 0.56-0.60 by Boroson ct al. 1996; (4) from the 
Hubble FOS at <j> ~ 0.5 by Anderson ct al. 1994. 



Table 5 
Absorption features 



Edge 


A (A) 


A (A) 


r 


State 


Neutral O 


23.05 


23.0 ±0.5 


0.30 ±0.06 


Low 




23.05 


23.27+°!° 


0.28 ±0.04 


Short-on 




23.05 


23.ll° i 


0.22 ±0.02 


Main-on 


Neutral N 


30.77 


(fixed) 


< 0.08 


Main-on 


Line 


A (A) 


A (A) 


EW (eV) 


State 


Atomic O ls-2p 


23.47 (1) 


23.50+1^ 


0.57 ±0.30 


Main-on 



Note. — Errors are 90% confidence limits. Symbols: A Q = theoretical edge wavelength; A = measured wavelength; r = edge optical depth; 
EW = equivalent width. 

References. — (1) McLaughlin & Kirby 1998. 
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Table 6 

Elemental abundance measurements and differential emission measure (DEM) parameters 



State 


(U/O) 


(N/O) 


(Me/O) 


K 

(10 58 cm" 3 ) 




X 2 of Fit 
(1 DOF) 


Ionizing 
Spectrum 


(C/O) 


(N/O) 


(No/0)g 


7 


Low 


0.52 ± 0.12 


8.8 ± 1.0 


1.9 ± 0.7 


7.2 ± 3.7 


1.59 ± 0.17 


0.02 


P.L.+cut (main-on) 


Short-on 


0.56 ± 0.15 


7.7 ± 0.9 


2.5 ± 0.7 


331- 


2.05 ± 0.19 


4.96 


P.L.+cut (main-on) 


Low 


0.30 ± 0.07 


5.9 ± 0.7 


2.5 ± 1.0 


10+| 


1.78 ± 0.17 


1.25 


20 keV brems. 


Short-on 


0.29 ± 0.08 


4.6 ± 0.6 


3.5 ± 1.0 


43+ 17 


2.28 ± 0.19 


8.72 


20 keV brems. 



Note. — We show the statistical errors at 66% confidence. The element abundance ratios (by number of atoms) are normalized to the solar 
values compiled by Wilms, Allen, & McCray (2000), which are (C/0)q = 0.47, (N/O) = O.lf , and (Ne/O) = 0.14. The K and 7 parameters 
define the DEM. 



